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ABSTRACT 

C 2 H is a representative hydrocarbon that is abundant and ubiquitous in the interstellar medium (ISM). 

To study its chemical properties, we present Submillimeter Array (SMA) observations of the C 2 H 
N= 3-2 and HC3N J=30-29 transitions and the 1.1 nun continuum emission toward four OB cluster¬ 
forming regions, AFGL490, ONI, W33Main, and G10.6-0.4, which cover a bolometric luminosity 
range of ~10 3 10 6 L@. We found that on large scales, the C 2 H emission traces the dense molecular 
envelope. However, for all observed sources, the peaks of C 2 H emission are offset by several times times 
10 4 AU from the peaks of 1.1 mm continuum emission, where the most luminous stars are located. By 
comparing the distribution and profiles of C 2 H hyperfine lines and the 1.1 mm continuum emission, 
we find that the C 2 H column density (and abundance) around the 1.1mm continuum peaks is lower 
than those in the ambient gas envelope. Chemical models suggest that C 2 H might be transformed to 
other species owing to increased temperature and density; thus, its reduced abundance could be the 
signpost of the heated molecular gas in the ^10 4 AU vicinity around the embedded high-mass stars. 

Our results support such theoretical prediction for centrally embedded ~10 3 - 10 6 Lq OB star-forming 
cores, while future higher-resolution observations are required to examine the C 2 H transformation 
around the localized sites of high-mass star formation. 

Keywords: astrochemistry; ISM: abundances; molecular processes; stars: early-type; stars: formation; 
stars: individual (G10.6-0.4, ON 1, W33 and AFGL490) 


1. INTRODUCTION 


Observations have shown that high-mass stars form 
in massive mole cular c lumps on scales of 0.1- 
1 


2009; 

Zhang et al. 2009; Wang et al. 2012t Liu et al. 

2011: 

Liu et al. 

2012allbl: Galvan-Madrid et al. 20131; 


Lu et all |2014). The strong pro- 
tostcllar feedback in radiative pressure and mechani¬ 
cal processes (e.g. outflows, stellar wind) starts to 
disturb the ambient molecular gas_whil e the mas¬ 

sive s tars are still deeply embedded dZinnecker fe Yorkel 
l2007fl . How the contracting clumps dynamically 
evolve under the influence of the (proto)ste llar feed¬ 
back remains uncer tain (e.g. iBeuther et all 120071 : 
iKrumholz et, al.ll2014fl . Extensive modelings have sug¬ 
geste d that this feedback may leave chemical footprints 
(e.g. lLahuis fe van Dishoeck l2000HChapman et alJl2009t 

iHerbst fe van Dishoeckll2009 : lTan et ah I2014H . High spa- 
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tial resolution observations of molecular lines from fast- 
evolving (10 4 - 10 6 yr) clouds are therefore helpful to di¬ 
agnose the history of cloud contraction and massive star 
formation even in the embedded phases. 

Ethynyl radical (C 2 H) is the simplest hydrocarbon 
molecule with the C=C bond, which is closely related to 
the formation and trans formati on of l ong carbon chains 
and other hydrocarbons dPilleri et al.ll2013h . C 2 H and its 
predecessor C 2 H 2 are crucial intermediates in the chem¬ 
istry of aromatic rings (e. g., CeHg) and polycyclic aro¬ 
matic hydrocarbon fPAH. lTielensll2013ll . However, be¬ 
cause of its symmetry, C 2 H 2 does not have a permanent 
dipole moment, so C 2 H 2 lacks of rotational transitions, 
which makes C 2 H 2 observationally challenging. On the 
other hand, C 2 H has many strong emission lines accessi¬ 
ble in millimeter bands. In addition, the hyperfine struc¬ 
ture lines (HFS) of C 2 H can provide estimates of opti¬ 
cal depth and magnetic field by measu ring the Zeem an 
effect . Its dipole moment (0.77 Debye, iWilson fe Greenl 
Il977f) is about seven times higher than that of CO, and it 
has a critical density of ~ 10 5 cm -3 (|Sakai et al.ll2010a l. 

C 2 H was first detected in Galactic star-forming re¬ 
gions bv lTucker et al.l (11974f) . who derived a comparable 
C 2 H abundance with those of HCN and HCO + . Later 
surveys found t hat CHI is widespread over the inner 
Galactic plane (iLisztl 119951 1. and it mainly arises from 
relatively dense g as (i.e. n(H 2 ) ~ 10 4 — 10 5 cm -3 , 

IWatt et all 11988 1. C 2 H has been detected in all 
evolutionary stages of high-mass star formation, e.g. 
infrared dark cloud (IRDC), prestellar cores, high- 
mass protostel lar objects (HMPOs) and hot molecular 
cores (HMCs) (iWootten et al.lfTOBOl : IBeuther et alJl2507t 
iLi et al.l I2012H . Moreover, extragalactic studies show 
that C 2 H is also quite strong, after those popular (dense) 
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gas tracers HCN, HNC, HCO + , CS, etc., so it can be 
used in multi-molecule diagnosis on galactic environ¬ 
ment (e.gJMeier fc Turneill2005ll2012HMeier et al.l[2015t 
LTiang et, a.1.1 1201 lb llmanishi fc Nakanishill2014fl . 

Single-dish and interferometric observations toward 
HMPO and ultracompact Hii (UCHii) regions found 
that C 2 H abundance is decreasing toward th e central hot 
cores dBeuther et al.ll2008l : iPilleri et al.ll2013h . which sug¬ 
gests that C 2 H can trace cold m olecu l ar gas in the early 
phase of star formation (see a lso iGerin et al • 20ll Ob¬ 
servations on prestellar cores (iPadovani et al.ll2009H also 
found similar trends of C 2 H abundance and implied that 
the sizes of C 2 H deplet ion holes are a few thousand AU. 
iMookeriea et all (|2012f) compared the distributions and 
abundances of C 2 H and C 3 in DR 21 (OH) and suggested 
the chemical (gas-phase) timescale for these species to 
be ~ 0.7 — 3 Myr. The aforementioned results, however, 
are still limited by spatial resolution and biased samples. 
A systematic survey will help clarify the environmental 
effects and will help us understand the formation and 
evolution of C 2 H in massive star-forming regions. 

In this paper, we present the Submillimeter Array 
(SMAjU observations of the C 2 H N = 3 — 2 transi¬ 
tions toward four massive star-forming regions, G10.6- 
0.4, W33 Main, ONI, and AFGL490. These samples 
were selected from the previous single-dish mapping sur¬ 
vey of C 2 H 1-0 (87.317 GHz), HNC 1-0 (90.664 GHz), 
and HC 3 N 10-9_[90.979 GHz) in 27 lrigh-mass star form¬ 
ing regions dLi et al.l [ 2012D , using the Purple Mountain 
Observatory 14m Telescope (PMO-14m; HPBW~55"). 
The accurate parallax distances have been measured ex¬ 
cept for AFGL 490, whose distance has an uncertainty of 
~50%. Details of our observations and data reduction 
are given in Section [ 2 ] We present the observational re¬ 
sults in Sections I3TT1 The analysis of the C 2 H optical 
depth and abundance profiles are provided in Section [3. 2 1 
and The implications of our results are discussed in 
Section [I] 

2. OBSERVATIONS AND DATA REDUCTION 

The SMA observations were carried out on 2012 Au¬ 
gust 4, with six available antennae in the subcompact 
configuration. The phase referencing and pointing cen¬ 
ters coincide with the wa ter m aser positions in the se¬ 
lected targets Isee lLi et al.ll2012t ). The weather condition 
was moderate (t 22 5 qhz ~ 0.25), yielding a typical system 
temperature T sys of ~300 K. The passband calibrator is 
3C 279, and the flux calibrators are Titan/Callisto. The 
absolute fluxes are accurate to within 15%. 

The observations covered a frequency range of about 
259-263 GHz in the lower sideband and 271 -275 GHz 
in the upper sideband. Our observations covered the 
C 2 H N= 3—^2, J= 7/2—>5/2, F=4-^3, i/ rest =262.00426 
GHz (blended with the F=3—>2 line), and J= 5/2—>-3/2, 
F=3—>2, i/ res t =262.06499 GHz (blended with the 
F=4—»3 line). These are the two main groups of the 
C 2 H N = 3 — 2 hyperfine structures. We simultane¬ 
ously obtained HC 3 N J = 30 — 29 (E u ~ 2 03 K) . which 
is an excellent tracer of HMCs (Chapman e t al.l 120091 

1 Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 


iMeier fe Turnerll2012D . A uniform spectral channel spac¬ 
ing was set to be 0.8125 MHz (~ 0.93 kms -1 ) across 
the entire passband. The SMA data were calibrated us¬ 
ing the MIR IDL software packagcQ We use the line-free 
channels in both the lower and upper sidebands to gen¬ 
erate the continuum channels, which are then imaged 
jointly to produce the 1.1 mm continuum image. The 
continuum emission is subtracted before performing the 
spectral li ne imaging. The imag ing is carried out using 
the CASA dMcMullin et al.l I2007T) package, with uniform 
weighting and hogborn CLEANing algorithm. We do not 
perform primary beam corrections since the angular sizes 
of the emission are small compared with the SMA pri¬ 
mary beam (~1'). In addition, the primary beam effects 
are factored out when we analyze the intensity ratios (see 
section I5T51) . Table (T] lists the observational properties of 
these sources, including the integration time, and the 
bolometric luminosities of these sources. 


3. RESULTS 

3.1. The C 2 H N=3-2 brightness distribution 

Figure [l] shows the integrated intensity maps of C 2 H 
N= 3-2 (E u ~ 25 K), which are overlaid with the 1.1 mm 
continuum emission (white contours) and HC 3 N <7=30- 
29 integrated intensity (cyan contours). The 1.1mm 
continuum emission may be contaminated by both the 
dust thermal emission and the free-free emission from 
the UC Hii regions. The liigh-excitation HC 3 N <7=30-29 
line (E u ~ 203 K) helps locate the HMCs. We resolved 
the extended C 2 H emission surrounding the continuum 
peaks (Figure [TJ, while the previous single-dish observa¬ 
tions of the lower excited C 2 H N=l-0 line ( E u ~ 4 K) 
traced th e diffu se emission from a more extended area 
dLi et al.l 120121) . The size and geometry of the warm 
molecular gas traced by C 2 H N = 3 — 2 are consistent 
with those traced by other molec ular lines (e.g. G10.6 - 
0.4: ILiu et al.llJOTfl W33Main: IHaschick fc HolU983lh 
suggesting that the resolved C 2 H N = 3 — 2 brightness 
distributions are likely not strongly biased by missing 
flux (Section [4| . The peaks of C 2 H N= 3-2 emission and 
continuum have systematic offsets . We do not r esolve 
ring -like C 2 H distributi on found in IBeuther et al.l (2008) 
and IGuelin et al.l (I1999H . likely because of the projection 
(or inclination) of the sample and the insufficient angular 
resolution of the observations. 

Figure [2] shows the intensity-weighted velocity maps. 
The C 2 H 3—2 emission shows velocity gradients, while 
its centroid velocities are comparable with the systemic 
velocities. More information about the individual targets 
is provided as follows. The measurements of the rms 
noise of the maps and the fluxes of C 2 H 3-2 are listed in 
Table [2] We achieve ~ 6 " angular resolution in the final 
images. 


( Ho & Haschick 

19861: Keto et al. 

119871 

1988 

Sollins et al.il 2005 

Liu et al.l 12010: Liu et al.l 

2010 


Liu et al] 1201 111 is a well-studied OB cluster-forming 


region ( Ho fc H aschickl 119811: ILiu et all I2013IL which 
is deeply embedded in the convergence of several ~5 
pc scale dense gas filaments. Adopting the parallax 


2 The MIR Cookbook by Chunhua Qi can be found at 
http://cfa-www.harvard.edu/~ cqi/mircook.html 
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Table 1 

General observational properties of the sample 


(1) 

Source 

(2) (3) 

R.A. (J2000) Dec 
(h m s) (° ' ") 

(4) 

Distance 

(kpc) 

(5) 

Luminosity 

(Lq) 

(6) 

Int.time 

(minutes) 

Gio.6-0.4 

18:10:28.70 

-19:55:48.7 

4.95 

7xl0 5 

120 

W33 Main 

18:14:13.67 

-17:55:25.2 

2.40 

5.4x10 s 

107 

ONI 

20:10:09.14 

31:31:37.4 

2.35 

2.7x 10 4 

120 

AFGL 490 

03:27:38.80 

58:47:00.0 

1.0 

2xl0 3 

100 


Note. — (1) Source names; (2) and (3) coordinates of the p hase centers; 
(4) Di stances of th e sources. Parallax d istances of G10.6 -0.4 from Sanna et al.l 
(120131) , W33 from llmmer et all (120131) , and ONI from IXu et al. 120130 : (5) 
Bolometric luminosities, converted from previous values using the latest par¬ 
allax distances, if available; (6) On-source integration time. 


Table 2 

Measurements 


(1) 

Source 

(2) 

Beam 

n 

(3) 

RMScont 
(Jy beam -1 ) 

(4) 

RMSc 2 h 

(Jy beam - 

(5) 

RMShc 3 n 
_ 1 kms —4 ) 

(6) 

■Sc 2 h 

(Jy kills -1 ) 

(7) 

Av a 

(km s —x ) 

Gio.6-0.4 

7.3x6.0 

0.20 

8.0 

2.0 

590 

[-10.7,2.8] 

W33 Main 

7.2x5.9 

0.12 

3.4 

2.6 

212 

[27.0, 40.5] 

ONI 

6.5x5.6 

0.04 

0.97 

1.3 

75.4 

[7.4, 14.6] 

AFGL 490 

8.9x5.0 

0.02 

0.83 


41.0 

[-17.1,-11.7] 


Note. — (1) Source names; (2) Synthesized beam (uniform-weighted) of the images; (3) 
1 cr RMS of the continuum images; (4) 1 a RMS of the C2H N=3-2 moment zero maps; (5) 
1 a RMS of the HC3N J= 30-29 moment zero maps; (6) Integrated flux of C2H N= 3-2; (7) 
Velocity ranges of the C2H 3-2 (7/2-5/2, 4-3) Main line. Note that the rms values listed here 
are higher than the expected thermal noise in the data since they are limited by the dynamic 
range in the images. 


distance of 4.95 kpc dSanna et al.l 120131 ). its bolometric 
luminosity is about 7 x 10 5 Lq (note that the luminosi¬ 
ties for all the sources in this paper are recalculated 
with their parallax distances, if available). 

On the ~0.5 pc scale, the C 2 H N = 3 — 2 emission 
(Figure 0 shows an obvious velocity gradient, and it 
likely traces the edge-on spinning-up molecular accretion 
flow resolved by the previous observations of highly ex¬ 
cited NH 3 lines dHo fc Haschickl H986t iKeto et all 119871 : 
ILiu et ahl 12010ft . The channel maps of C 2 H N = 
3 — 2 (Figure 0 demonstrate the known southeast- 
north west velocity gradient in this fl attened accretion 
flow (IKeto et al.Hl987t ILiu et all 12010ft . The integrated 
intensity of C 2 H peaks at both blueshifted and redshifted 
lobes, which locate to the SE and NW. It shows weaker 
emission at the center. We also notice that the C 2 H emis¬ 
sion is weaker in the west of the region, where the other 
two UCH 11 regions reside (red asterisks in Figure [I] and 
0 . 

W33Main — W33Main (G12.8-0.2) resides in_the OB 
clust er-forming region W33 Complex (Haschick fe Hoi 
11983ft . At a distance of 2.4 kpc (llmmer et al.l l2013fT T 
its bolometric luminosity is about 5.4 x 10 5 Lq. In our 
observations, the 1 . 1 mm continuum (~3 Jy) peaks at 
R.A.(2000) = 18 h 14 m 13.8 s , decl.(2000) = -17°55'43.4", 
which is assoc i ated with the infrared source IRS 3 
(|Dvck fc Simonl 11977ft and the centimeter peak W33 
Main-B2. The distribution of HC 3 N 30-29, however, 
is slightly different and is closer to anot her centimeter 
emission peak, B1 (|Haschick fe Hoi [19831. We also de¬ 
tect another 1.1 mm source associated with the IRS 2 and 


the UCH 11 region W33 Main-A in the north. It shows 
weaker emission in 1.1 mm continuum (~ 0.6 Jy) and 
in the HC 3 N 30-29 transition, comparing to the south 
component. 

As shown in Figure [D the C 2 H N = 3 — 2 emission is 
mainly distributed around W33 Main-A and W33 Main- 
B. The emission of C 2 H 3-2, HC 3 N 30-29 and 1.1mm 
continuum around W33Main-B are brighter and appear 
to be offset from each other, while in W33 Main-A their 
peak positions agree better with each other. The envi¬ 
ronment of W33 Main is quite complicated since there 
are multiple (UC)Hii regions, and these (UC)Hii re¬ 
gions A, Bl, and B2 seem to be embedded in the gas 
envelope traced by C 2 H 3-2. Although the projection 
effect cannot be ruled out, the differences between these 
two regions (W33 Main-A and B) could be explained by 
their different physical conditions and star-forming ac¬ 
tivities. Comparing with W33 Main-A, W33Main-B is 
much more massive and brighter in infrared, millimeter, 
and centime ter continuum, as wel l as consisting of more 
Hii regions ((Haschick fc Holfl983ft . As revealed by the 
continuum and HC 3 N on ~ 0.1 pc scales in Figure 0 
the different behaviors of Bl and B2 also suggest that 
these two regions have different stages of SF process. 
The velocity field in W33 Main-B is complicated (Figure 
0 , which sug gest s infl uence from the multiple molecular 
outflows (llmmer et al.ll2014 h 

ONI — The UCH 11 region Onsalal (ON 1) is a cluster- 
forming region with multiple sta r-forming signposts 
(|Kumar et al.1120031: ISu et all 12009ft . At a parallax dis¬ 
tance of 2.35 kpc (|Xu et al.ll2013ft its luminosity is about 
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Figure 1. Grayscale and black contours are the C 2 H N = 3 — 2 (J =7/2-5/2, F = 4 — 3) uniformly weighted moment zero (integrated 
intensity) maps of G10.6-0.4, W33, ONI and AFGL490, respectively. C 2 H intensity levels (black contours) start from the 5cr level and 
continue in 3 a steps. The white contours are the 1.1 mm continuum images (30%, 50%, 70%, and 90% the peak emission). The continuum 
peaks are 8.0, 4.1, 1.8, and 2.0 Jy beam -1 , respectively. The cyan contours are HC 3 N J — 30 — 29 uniformly weighted moment zero maps. 
In G10.6-0.4 the HC 3 N 30-29 contours start from the 5 cr level and continue in 10 a steps. In W33Main and ON 1 the HC 3 N 30-29 contours 
start from the 4<r level and continue in 3cr steps. Th e 1 cr of the continuum, C 2 H 3-2 and HC 3 N 30—29 are listed i n Table [2] The re d 
asterisks denote the positions of the (UC)Hn regions dHo Sz Haschick 1986; Hasc hlck Sz Hdll983t [Kumar et al.N2004l ; ISimon et alj [ l9 83). 
The yellow crosses denote the position of the H 2 O masers. The SMA subcompact synthesis beams of the C 2 H images are shown in the 
bottom left corners. A colorful figure is available in the online version. 


2.7 x 10 4 Lq IlMacLeod et al.llT998l l. Figure Q] shows that 
the 1.1 mm continuum and the HC 3 N J = 30 — 29 emis¬ 
sion both agree well with the center of the UC Hu and the 
H 2 O maser, while the C 2 H N=3-2 emission peaks show 
offsets in the southwest. Figure [2] shows that C 2 H has 
a velocity gradient along the northeast-southwest direc¬ 
tion, and the velocity at the emission peak is c onsistent 
with the systemic velocity vlsr =12 kms -1 feu et ali 
I2009T) . Such distribution along the NE-SW direction 


is in agreem ent with the H 13 CO + distribution reported 
by iKumar et al.l (12004( 1 , which is interpreted as one of 
the multiple molecular outflows in ON 1, while high- 
res olution obse rvations toward the innermost region by 
ISu et alJ (I2009H favor a scenario of the expansion of the 
Hii region. Other than the main SW component, Fig¬ 
ures [T] and [ 2 ] also show that the C 2 H emission appears to 
slightly extend to the northwest, in a direction perpen¬ 
dicular to the proposed bipolar outflow (or expansion 
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Figure 2. Moment one (intensity-weighted average velocity) maps of C 2 H N = 3 — 2 (J =7/2—5/2, F=4-3) of G10.6-0.4, W33, ON 1 and 
AFGL490. The black contours are the C 2 H 3-2 moment zero, and the white contours are the 1.1 mm continuum, as shown in Figure 1X1 
The asterisks and crosses denote the (UC) Hu regions and the H 2 0 masers. The SMA subcompact synthesis beams (~ 6") are shown in 
the bottom left corners. A colorful figure is available in the online version. 


motion). This component might be related to the other 
outflow com ponent traced by SiO and CO as reported by 
iKumar et al.1 (|2004f) . These authors derived that the cen¬ 
tral enclosed mass of ON 1 is about 300-400 Mq, which 
would be sufficient to support a rotating structure, so 
with current data we are not yet able to rule out the 
possibility of a flattened rotating accretion flow confin¬ 
ing the UC Hu region. 

AFGL490 — AFGL490 is a deeply embedded 8-10 M 0 
young star. Previous CO and 13 CO observations sug¬ 
gest a warm (Tki n >100 K) gas component within the 
central ~ 6000 AU radius, surrounded by a cooler gas 


envelope dMitchell et al.lfl995t ). iSchrever et al.1 (120021) re¬ 
ported that the size of the outer envelope is about 22,000 
AU x 6000 AU. Molecular outflows have also been in¬ 
ferred from the detection of high-velocity CO emission 
(|Mitchell et al.l[T995h . 

In AFGL490 C 2 H 3-2 is significantly weaker than that 
in other sources, and HC 3 N 30-29 is not detected. The 
brightness distribution of C 2 H N = 3 — 2 mainly fol¬ 
lows the flattened dense gas envelop e traced by the CS 
.7 = 3 — 2 transition (Figure fTl ISchrever et al.ll 2002 D and 
peaks at northwest of the 1.1 mm continuum source. The 
velocity gradient of C 2 H N= 3-2 in Figure [2] is not as 
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Figure 3. Examples of C 2 H N = 3 — 2 channel maps of G10.6-0.4, 
from Vlsr — "0-^ kms -1 to -3.5 kms -1 . Contours start from 5cr 
and increase in 3 cr steps (dashed contours are negative). The SMA 
subcompact synthesis beams are shown in the bottom left corners. 
A solid square in the first channel map denotes the position where 
we extract the spectrum in the bottom panel. The spectrum shows 
two velocity components of C 2 H 3-2 emission in that position. 
Numbers “1” and “2” denote the main hyperfine groups of the two 
velocity components, and the red dashed line is the fitting to the 
spectrum. The fitting residual and corresponding RMS are shown 
below the spectrum. A colorful figure is available in the online 
version. 


clear as in G10.6-0.4 or ON 1, but it still shows that the 
C 2 H N= 3-2 emission consistently traces the northwest- 
southeast velocity gradient in this flattened dense gas en¬ 
velope, as well as the molecula r ga s wind perpendicular 
to the flattened envelope dSchrever et all[20021 . 

3.2. Optical depth of C 2 H N = 3 — 2 

The electron-nucleus interaction of C 2 H results in six 
splitting hyperfine line components. Optical depth can 
be obtained by fitting the hyperfine structures. We ex¬ 
tract the C 2 H N = 3 — 2 spectra in every ~3" x3" pixels 
and fit the optical depth using the CLASS package^ based 


on the the oretical intens ity ratios provided in the CDMS 
database^] (iMiiller et al Jl2005h . LTE relative line intensi¬ 
ties were assumed. The bottom panel of Figure [3] shows 
an example spectrum of G10.6-0.4 extracted from a po¬ 
sition offset from the core center (indicated by the solid 
square in the first channel map of Fig. ©. The spectrum 
shows two velocity components (line separation AV ~ 5 
kms -1 , labeled T’ and ‘2’), which would cause confusion 
in the fitting of hyperfine structures (the velocity sepa¬ 
ration of the two hyperfine lines of each group is about 
2.54 kms -1 ). We use Gaussian profiles to decompose 
the two velocity components, and then we derive their 
optical depths separately. 

We found that C 2 H is optically thin (r ~ 0.1) in most 
regions, which is comparable with the optical dep th rf 
C 2 H found in the UCH 11 region MonR2 dPilleri et al.l 
120m Using a single-dish telescope iLi et al.l ([2012D also 
found that the optical depths of C 2 H are small to moder¬ 
ate in most massive star-forming regions. The estimates 
of the C 2 H column density ./Vc 2 h and abundance A'c 2 h 
are given in section! 



Figure 4. Normalized intensity profiles of different tracers (red 
solid lines-. C 2 H 3-2, blue dashed lines: continuum; magenta dotted 
lines: HC 3 N 30-29) of the four sources, which are extracted and 
measured along a strip starting from the continuum peaks to the 
C 2 H peaks. Horizontal axis is the distance from the continuum 
peaks. A colorful figure is available in the online version. 


3.3. Intensity and column density profiles 

The HC 3 N J = 30 — 29 emission (E u ~ 203 K) provides 
a good probe of the dense and warm gas concentration. 
We did not detect HC 3 N 30-29 in AFGL490. For the 
rest of the sources, the normalized 1.1 mm continuum 
and the radial distributions of C 2 H 3-2 and HC 3 N 30- 
29 are presented in Figure [H We use one strip from the 
continuum peak to the C 2 H peak, to extract the intensity 
profile for each source, so as to show the relative offset 
scale of different tracers. 

The C 2 H 3-2 intensity increases outward until sev¬ 
eral times 10 4 AU separation from the 1.1 mm emission 
peak and then drops steeply. In G10.6-0.4 and ON 1, 


3 http://www.iram.fr/IRAMFR/GILDAS 


4 http://www.astro.uni-koeln.de/cdms 
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HC 3 N 30-29 appears to be correlated with the 1.1 mm 
continuum emission and is confined in the compact re¬ 
gions (see also Figure [TJ. In W33Main, HC 3 N 30-29 
peaks in between the peaks of the 1.1 mm continuum and 
the C 2 H 3-2 emission, and HC 3 N is very close to one of 
the Hii region B1 (Figure[U). Nevertheless, the C 2 H 3-2 
emission peak is spatially offset from the 1.1 mm contin¬ 
uum and the HC 3 N 30-29 emission peak in W33 Main. 

To quantify the column density and abundance pro¬ 
files of C 2 H, we extract the C 2 H 3-2 line intensity and 
the 1 . 1 mm continuum flux in several regions with di¬ 
ameter ~ Z" from the continuum peak and through the 
C 2 H 3-2 line emission peak, until the S/N drops to < 2 a 
for either the 1.1mm continuum or the C 2 H 3-2 emis¬ 
sion. Column densi ties of H? are derived from 1.1mm 
continuum, following iBeuther et al.l (120021(2003 1 and as¬ 
suming an_emissivit^_^_= 1.5 a nd fti. 2mm = 0.97. Fol¬ 
lowing iHunter fc Broganl (12006H , the column density of 
C 2 H, A(C 2 H), is derived using 


AT(C 2 H) 


Zirkv 2 Q t 
hc 3 A ui g u 


exp( 


E u 

kT ex 


) ■ -S'c 2 h 


(1) 


where A u i, Qt.,gu are the Einstein A coefficient, the par¬ 
tition function, and the upper-state degeneracy quoted 
from the CDMS, respectively. 5 c 2 h is the line inten¬ 
sity of C 2 H in units of Kkms -1 . Identical T ex = 20 
K and Tdust = 50 K are adopted in the calculations for 
all the sources. Such an assumption on Tdust is consis¬ 
tent with previous values adopted for these sources (e.g. , 
Liu et al.ll2010l: lHaschick fc Holll983l: iKumar et alJl2004 

Mueller et al.ll2002t iSchrever et al .1120001 1. where the T ex 
are found to be all close to 20 K in the hyperfine struc¬ 
ture fittings of C 2 H (Section 13.21) , similar to the value 
derived by iDe Beck et al.1 (|2012l h We find that the de¬ 
rived iV(C 2 H) values are a few times 10 15 cm ~ 2 for all 
targets. 

There are several sources of errors in the estimates of 
the C 2 H column density and abundance relative to the 
H 2 molecule. The variation of T ex between 15 and 25 K 
can introduce up to 0.3 dex error in log 1V(C 2 H), so the 
typical error in log 1V(C 2 H) is about 0.4 dex considering 
the error on Sc 2 h of about 50%. The variation of Td us t 
between 40 and 60 K would also result in about 0.1 dex 
error in log _/V(H 2 ). Although adopting different values of 
0 (e.g., 0 = 1, 1.5, 2) will introduce about 0.3 dex uncer¬ 
tainty in log A(H 2 ), we expect that 0 is probably similar 
in these sources. Uncertainties on the absolute fluxes of 
continuum and C 2 H line are factored out as A(C 2 H) is 
based on their intensity ratios. Thus, the typical error in 
log X (C 2 H) is about 0.5 dex. However, we caution that 
the derived C 2 H column density and abundance profiles 
cannot be very precise at this moment. This is mainly be¬ 
cause ( 1 ) part of the 1.1 mm continuum may be contam¬ 
inated by free-free emission, which cannot be accurately 
subtracted yet, and ( 2 ) the missing fluxes in interfero¬ 
metric observations can lead to uneven underestimates 
of the continuum and line intensities. Thus, the results 
could be possibly affected by another factor of two. Nev¬ 
ertheless, as we used the SMA subcompact configuration 
with about 6" resolution, thos e exten d ed em issions on a 
scale of less than 1' (~ 1 pc; iLi et al.l I2012T) should be 
recovered effectively, so we suppose that the amount of 
diffuse emissions should be limited and the results are 




Figure 5. Top: C 2 H column density radial profiles in logarithm 
scale of the four sources. Horizontal axis is the distance from the 
continuum peaks in units of AU. The typical error on the AT(C 2 H) 
is about 0.4 dex. The horizontal lines represent the iV(C2H) of 
other types of star forming regions (Gerner ct al. 2014; Sakai et al. 
2010b). Bottom: C 2 H abundance profiles. The points with ar¬ 
rows represent lower limits of the C 2 H abundances. The error bar 
in the bottom right corner denotes the typical error on the C 2 H 
abundances. On the top of the figure, the horizonta l dashed lines 
represent the C 2 H abundances from oth er works d Gerner et al.l 
120141 : [Lee et alJ 1 19981 : [Lucas fc LisztifTOOOl) . The diamonds on the 
uppe r left are the C 2 H abundance profiles of the prestellar cores 
dPadovani et al.l120091) . Colorful figures are available in the online 
version. 


not affected significantly. 

Figure [5] shows the A r (C 2 H) and A'(C 2 H) as a function 
of the separation (physical scale) from the 1.1 mm contin¬ 
uum peaks. For comparison, we also plot some A(C 2 H) 
and A(C 2 H) from other works in Figured] On the top 
panel, the horizontal reference lines denote the median 
A(C 2 H) of differ ent types of star formation regions ob¬ 
tained bv I Gerner et alJ (|20l4 dashed), the i V(C 2 H) of a 
low-m ass star-forming region L1527 (dotted, iSakai e t alJ 
l2010bf) . and TMC 1 (dot dashed, lOhishi et al.l 11995b 
Three of our sources have similar C 2 H column densities 
as the HM PO, but sligh t ly low er than the N(C 2 H) of the 
UCHu in iGerner et all (2014 b Note that the reference 
lines trace very different spatial scales as they represent 
single-dish observations. The _/V(C 2 H) of AFGL490 is 
lower than other sources, which might be due to the fact 
that the source is at a latter evolutionary stage when 
many molecules are destroyed by the ionizing radiation. 
For G10.6-0.4 and W33Main the A(C 2 H) in the centers 
are about a factor of two lower than their peak values, 
and overall the decreasing trend of _/V(C 2 H) in the core 
centers of our targets is not very obvious. 

The bottom panel of Figure [5] shows the X(C 2 H) as 
a function of radius. We found that the A'(C 2 H) at 
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Figure 6. Top: C 2 H column density, 7V(C2H), as a function of 
bolometric luminosity. The symbols and reference lines are the 
same as in Figure [51 The typical error on the 7V(C2H) is about 0.4 
dex. Bottom: C 2 H abundance X(C 2 H) as a function of bolometric 
luminosity. The diamonds on t he upper left are the C 2 H abundance 
of the prestellar core s L149 8 {Padovani ct al. 2009]), and its Lbol 
is from Shirley ct al. (2005). Colorful figures are available in the 
online version. 


larger separation are close to the values in the prestel- 
lar co res (~10~ 8 , denoted as diamonds; iPadovani et al.l 
120091 1 but lower than those of other types of clouds (ref¬ 
erence lines). X(C 2 H) systematically decreases toward 
the 1.1 mm continuum peaks by about one order of mag¬ 
nitude, and the four sources have similar slopes of the 
C 2 H abundance radial profiles. Comparing the plots of 
./V(C 2 H) and A'(C 2 H) together we find that our derived 
A'(C 2 H) are lower than that of other work, and the in¬ 
consistency between our data and other work is mainly 
for the abundance plot. 

Figure p plots iV(C 2 H) and X(C 2 H) as a function of 
bolometric luminosities of the four sources. The top 
panel shows a trend that iV(C 2 H) increases in more lu¬ 
minous targets, and the more luminous sources G10.6- 
0.4 and W33Main have the largest 1V(C2H), despite the 
small sample and the uncertainty on iV(C 2 H). While the 
luminosities of the sources differ by nearly three orders 
of magnitude, their mean 1 V(C 2 H) differ only by about 
a factor of four, (see Section [4] for more discussion). On 
the other hand, the bottom panel of Figure [G] shows that 
AT(C 2 H) of our sample seem to be lower than the less 
luminous source L1498, which is a prestellar core. We 
again note that the uncertainty of A(C 2 H) is large, al¬ 
though there is a plausible trend in the bottom panel of 
Figure [5] and [G] 

4. DISCUSSION AND CONCLUSIONS 

We used the SMA to observe the C 2 H N= 3-2 emission 
lines (~262 GHz) toward four high-mass star-forming re¬ 
gions G10.6-0.4, W33Main, ON 1 and AFGL490 with a 
spatial resolution better than 0.15 pc scale. Our achieved 


resolution is comparably smaller than the typical pro¬ 
jected diameters of the HMCs, which is the key to exam¬ 
ining the C 2 H chemical evolution under the (proto)stellar 
illumination. Our results show that the extended C 2 H 
emission traces the dense gas envelopes in these high- 
mass sources, and its column density V(C 2 H) is of a few 
10 15 cm -2 . 

Our observations show that, on ~ 0.15 pc scales, the 
C 2 H emission peaks are systematically offset from the 
high-mass star-forming hot cores. This is consistent with 
the previo us S MA observation of the source IRAS 18089- 
1732 dBeuther et al.l[2008l ), as well as their model, which 
interprets such a scenario as the transformation of C 2 H 
into other molecules in the dense and hot environment 
of the core center. The profiles in Figure 0] support such 
a proposed picture. Our findings of high column den¬ 
sity C 2 H are also in agr eem ent with the single-dish sur¬ 
vey by iGerner et al.l (|2014f) , in more luminous sources 
(Figure [6]). This might be a natural result that more 
luminous sources have higher gas masses and a stronger 
PDR (photo-dominated region) e ffect. The age of the 
specific HMPO model presented in lBeuther et al.l (120081 ! 
was 5xl0 4 yr, but since the ages of our observed sources 
are uncertain, and the abundance of C 2 H cannot be de¬ 
rived very accurately at this moment, a direct compar¬ 
ison between the models and the observations remains 
difficult. 

In Figures [5] and [6] we compare our derived X(C 2 H) 
with that from the previous observations toward a sam¬ 
ple of high-mass star-forming regions at di fferent ev olu- 
tionary stages (HMPO, HMC, and UC Hu. IGerner et al.1 
[Ml .' All sources Except for AFGL490 show simi¬ 
lar V(C 2 H), which are close to the median Ai(C 2 H) of 
the HMPO, but slightly lower than that of t he U CH 11 . 
We note that in the work by IGerner et al.l 1 2014 1 the 
C 2 H data of the HMCs and UC Hii (which are treated 
as later stages) were not well explained by the model, 
which might imply a more complicated chemistry of C 2 H. 
X(C 2 H) of the four sources are also comp ared w ith that 
of other clouds, including that from IGerner et all ( 2014 ). 
the dark clouds TMC-1 and L134N, the prestellar cores, 
and the diffuse molecular clouds (Section EOl) . The abun¬ 
dances derived for our targets are significantly lower than 
those from other clouds, which is likely induced by the 
uncertainty of the interferometric observation. BuUnote 
that the X(C 2 H) of the prestellar cores (Pa dovani et al.1 
2009) are also lower than in those works, and their resolu¬ 
tions are comparable to the SMA (~ 0.15 pc) as their tar¬ 
gets are much closer, so we suspect that the low X(C 2 H) 
we find could also be a result of the difference between 
the physical resolution of the works. 

Figure [6] shows an increasing trend of _/V(C 2 H) in more 
luminous sources, while in AFGL 490 the mean X(C 2 H) 
is the lowest. Although AFGL 490 is still embedded in a 
molecular envel op, it is already in a transition stage to 
Herbig Be stars dSclirever et al.ll2006l! . so it is likely that 
the neutral molecular gas content of AFGL 490 has been 
reduced as a result of the ionized radiation. In fact, the 
reason for this trend in Figure [6] is likely that more lumi¬ 
nous sources tend to have more molecular gas, but as our 
images (Figure [Tj) show, while C 2 H is weaker in the cen¬ 
tral regions where it is possibly transformed, this trend 
show that the overall C 2 H mass can maintain a certain 
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amount in the gas envelop surrounding the cores. This is 
also consistent with the picture that C 2 H is a good tracer 
of PDR, in the sense that it trace s the cloud surface illu - 
minated by the stellar radiation (Meier fc Turnedl2005l h 
It is also possible to compare this trend with that of 
other molecular lines to analyze the corresponding chem¬ 
ical process or radiation effect on molecules. Owing to 
the uncertainties mentioned above that cannot be accu¬ 
rately accounted for at this moment, the quantitative 
analysis of A(C 2 H) is less conclusive, and it is unclear 
yet whether the abundance of C 2 H is significantly af¬ 
fected by the bolometric luminosity. Short-spacing data 
are needed to improve the accuracy of the abundance, 
and higher-resolution observations can help to discern 
the genuine spatial structure of C 2 H in different environ¬ 
ments, as well as to provide additional information on 
the < 10 4 AU scales where models of different parame¬ 
ters start to deviate fiBeuther et al.l 12008 1. 

4.1. Chemistry of hydrocarbons 

Observations toward the PDR region MonR2 found 
that C 2 H is extended compared to CH and c-C 3 H 2 , and 
the abundance of C 2 H is constant over a wide range of 
Go (UV intensity; iPilleri et al.112013 1. They demon¬ 
strated that the behavior of C 2 H is dominated by time- 
dependent effects and gas-grain chemistry in low-density 
molecular envelopes. This scenario agrees with our ob¬ 
servations, which also show evidence that C 2 H is stronger 
in the molecular envelopes. 

The formation of C 2 H is mainly determined by the 
formation and evolution of ionized hydrocarbons, which 
involves the association of C + and H 2 : C + —1 CH+ 
(n = 2, 3,4, 5). CH+ can be transformed into CH 4 , then 
C 2 H+ and C 2 H^". C 2 H is formed from the dissociated re- 


combination of C?H.t and CqH^ (see, e.g., 

Turner et al.1 

2000; Meier & Turned 2005); Beuther et al. 

200 S for de- 

tails): 


C‘2^2 e ^ C 2 H + H, 

( 2 ) 

C 2 H+ + e" -aC 2 H + 2H/H 2 , 

(3) 

Other pathways of producing C 2 H include the pho¬ 
todissociation of acetylene (C 2 H 2 ), and the reaction be- 

tween CH 2 and carbon atoms (Sakai et al. 20101J): 

C + CH 2 -+ C 2 H + H. 

(4) 


On the other hand, C 2 H is mainly destroyed by the 
reactions with oxygen atoms: 

C 2 H + 0—1 CO + CH. (5) 

It can also be transformed into a heavier carbon chain 
such as C 5 (|Mookeriea et al.ll2012H . Models addressing 
the variation of ga,s-phas e [C]/[0] abundance ratio (e.g., 
iTerzieva fe Herbstl Il998f) might help constrain the pro¬ 
cess of production/depletion of C 2 H in these sources. 

It would be also interesting to compare C 2 H with other 
simple carbon chain molecules such as I K AN, s i nce t hey 
share the same chemical precursor C 2 H;)" (iWootten et all 
ll98d ). Our observation covered HC 3 N J = 30 — 29 simul¬ 
taneously but its high excitation energy ( E u ~ 200 K) 
makes it difficult to compare its chemistry with C 2 H 
N = 3 — 2 (E u ~ 25 K) because they trace the dis¬ 
tinct physical conditions. Observations at 3 mm band 


would be able to cover both C 2 H 1-0 (87.4 GHz) and 
HC 3 N 10-9 (91 GHz, E u ~ 24 K) simultaneously, and 
such comparison is more suitable for studying the chem¬ 
istry of carbon chain molecules in cold gas. The compar¬ 
ison between C 2 H and other molecular lines and related 
chemistry will be further addressed in future work. 

4.2. Extragalactic perspective 

Our derived _/V(C 2 H) can also be compared with stud¬ 
ies on external galaxies. Obs ervation of 1 kp c scale 
molecular complex in M 51 by IWatanabe et all (I2014T ) 
found that _/V(C 2 H) is a few times 10 14 cm -2 (emis¬ 
sion size corrected, corresponding A(C 2 H) is about ~ 
4.5 — 5 x 10~ 9 ), which is about an order of magnitude 
lower than the average A(C 2 H) of the sources in our 
work. This might suggest that while C 2 H has higher 
column density near PDRs on ~ pc scales, and observa¬ 
tions on large scales (~ 1 kpc) would mainly collect C 2 H 
emission from those diffuse molecular gas._ 

In the high - resolution line surv eys bv lMeier fc Turned 
B f2012h : iMeier et all (l20T5h . C 2 H appears to be 
bright (comparable to CS in NGC 253) and well traces 
the circunmuclear PDR regions of the nearby galaxies, 
and even traces the molecular o utflow in the nucleus of 
Maffei 2 (jMeier fc Turned l2012fl . The hyperfine struc¬ 
tures of C 2 H can help estimate the optical depth, and 
the frequency of C 2 H 1-0 is very close to HCN 1-0 and 
HCO + 1-0, so they can be observed simultaneously, and 
their different critical densities and chemical properties 
can be c om pared to diagnose the environment of galax¬ 
ies. I.Iiang et aD ( 201 if have reported a sample of active 
galaxies detected in C 2 H with the IRAM 30m telescope, 
and with the aid of ALMA and other facilities, the role 
of C 2 H in studying extragalactic molecular gas will be 
further explored. 
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